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MORSE FUNCTIONS OF THE TWO SPHERE 


LIVIU I. NICOLAESCU 


Abstract. We count how many ’’different” Morse functions exist on the 2-sphere. There 
are several ways of declaring that two Morse functions / and g are ’’indistinguishable” 
but we concentrate only two natural equivalence relations: homological (when the regular 
sublevel sets / and g have identical Betti numbers), and geometric (when / is obtained 
from g via global, orientation preserving changes of coordinates on and R). The count 
of homological classes is reduced to a count of lattice paths confined to the first quadrant. 
The count of geometric classes is reduced to a count of certain labelled trees. We produce a 
two-parameter recurrence which can be encoded by a hrst order quasilinear pde. We solve 
this equation using the classical method of characteristics and we produce a closed form 
description of the exponential generating function of the numbers of geometric classes. 
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1. The main problem 

Suppose X is a smooth compact, oriented manifold without boundary. Following R. 
Thom, we say that / ; X ^ M is an excellent Morse function if all its critical points are 
nondegenerate and no two of them correspond to the same critical value. We denote by 
Mx the space of excellent Morse functions on X. To ease the presentation, in the sequel a 
Morse function will by default be excellent. 

For / G Mx we denote by p(/) the number of critical points of /. Given a Morse 
function / : X ^ M with n{f) = n we define a slicing of / to be an increasing sequence of 
real numbers 

— OO = ao < Ul < • • • < Un-l < CLn = oo 

such that for every i = 1, ■ ■ ■ ,n the interval (oj-i, a*) contains precisely one critical value 
of /. 

Two Morse functions f,g:X^'R. will be called geometrically equivalent if there ex¬ 
ists an orientation preserving diffeomorphism r : X —> X and an orientation preserving 
diffeomorphism : R ^ M such that 

g = £ o f o r. 

We denote by this equivalence relation. 

Two Morse functions /, </ : X ^ R will be called topologically equivalent if p(/) = ^{g) = 
n and there exists a slicing oi < • • • < On-i of /, a slicing hi < ■ ■ ■ < hn-i of g and 
orientation preserving diffeomorphisms 

: {/ < ai} ^ {g < h}, Vi = 1, • • • , n, = 6n = oo. 

Two Morse functions /, </ : X —> R will be called homologically equivalent if p(/) = y{g) = n 
and there exists a slicing ai < • • • < On-i of / and a slicing 6i < • • • < bn-i of g such that 
for every i = 1, • • • n the sublevel sets {/ < Uj} and {^f < hi} have the same Betti numbers. 
We denote by ~t and these equivalence relations. Note that 

/o /i /o fi fo /i- 
Set 

{/ V Mx; p(/) = n}, [Mx]* := Jdx/ * € {h,t,g}. 

Observe that we have natural projections 

[Mx]g ^ P^x]t P^x]h- 

Since the excellent Morse functions are stable we deduce that the geometric equivalence 
classes are open subsets of M[^. This shows that the quotient topology on [M^^]*, * G 
{g, t, h}, is discrete. 

If we think of a Morse function / G M[^ as defining a sort of “triangulation” on X with 
n simplices then we can expect that [M^j^ is finite^ and therefore the sets [M^jt^/i ought 
to be finite as well. Clearly, the sets M[^ contain information about the manifold X. For 
example, Morse theory implies that 

M^^0^n> j;6fc(X). 

k 

We can then ask how big are of the sets [M^^]*, * G {g,t,h} and what kind of additional 
structure do they posses. It is more realistic to start by answering this question for special 

haven’t worked out a formal argument showing that [Mx]" is hnite, but I believe this to be the case. 
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X. The simplest case is X = 5'^. We discuss briefly this case since it bares some similarities 
with the case X = we will discuss at length in this paper. 

The sublevel sets of a Morse function on are disjoint unions of closed intervals and we 

deduce that in this case the topological and homological classifications coincide. A Morse 
function on has an even number of critical points, but two of them are very special, 
namely the global minimum and maximum. 

Suppose the Morse function / : 5^ —> M has 2m+ 2 critical points Po,pi, - ■ ■ ,p 2 m+i, with 

Cj = f{Pi) < f{Pj) = Cj iff f < j. Set Li := {/ < Cj + e}. Then Lq and L 2 m are closed 

intervals. If we set Xj(/) := hQ{Li) then we see that a Morse function defines a sequence 

x{f) : {0,1, 2, • • • , 2m} ^ Z>o, i i—> Xi(/), 

satisfying 

xq = X2m = 1, Xi> 0, |xi+i - Xi\ = 1, Mi. 

Conversely, to any such sequence we can associate a Morse function with 2m + 2 critical 
points and we have 

/ 5 x{f) = x{g). 

We regard such a sequence as a walk of length 2m on the lattice Z C M with steps of size 
±1 which starts and ends at 1 and it is confined to the positive chamber Z>o. The number 
of such walks can be easily determined using Andre’s reflection principle, [Tni Ex. 14.8] and 
we obtain 



where Ck = i^k) denotes the k-ih. Catalan number. 

The number of geometric equivalence classes of Morse functions on can be determined 
using the calculus of snakes of Arnold, |2ll2j- We outline the main idea. 

Suppose / : 5^ —> M is Morse function with 2m + 2 critical points. Denote by p-i G 
the point where / achieves its global maximum. Then, starting at p-i label the criti¬ 
cal points PO) ■ • • )P 2 m, in counterclockwise order (see Figure P). Now remove the arcs 
[p_i,Po)) (P 2 m,P-i] (see Figure P). What’s left is what Arnold calls a A 2 m-snake. Their 
numbers are and the associated generating series are determined in HO]. 



Figure 1 . Cutting a Morse function down to a snake. 

Let us point out that the geometric classification of Morse functions on is different 
from the topological classification. For example the Morse functions depicted in Figure H 
are topologically equivalent, yet not geometrically so. In this figure the Morse function is 
the height function y, the vertices are the critical points and the numbers attached to them 
are the corresponding critical values. 
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Note that for the function on the left hand side the first and third critical points lie in 
the same component of the sublevel set {y < 4 + e} but this is not the case for the function 
on the right hand side. 
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Figure 2. Two Morse functions which are topologically equivalent but not 
geometrically equivalent. 


2. The main results 

In this paper we investigate the structure of the sets where * E {h,g}. Here is 

briefly its content. 

In Section 3, using the basics of Morse theory we analyze what kind of handle additions 
can occur as we cross the critical levels of a Morse function on S'^. This leads immediately to 
a bijection between the homology classes of such Morse functions and sets of lattice paths 
in The cardinalities of such sets are computed explicitly in Theorem 14.11 using the 

reflection technique of Gessel and Zeilberger, [Jj. More precisely, in Theorem 14.1 1 we prove 
that the number of homology classes of Morse functions on 5^ with 2n critical points is 

h2n — CnCn—1: 

where Cn denotes the n-th Catalan number. 

In Section 5 we describe a bijection between the set of geometric equivalence classes of 
Morse functions on 5^ and the set of Morse trees. These are labelled trees with vertices 
of degree one or three such that every node (i.e. vertex of degree three) has at least 
one neighbor with a bigger label and at least one neighbor with a smaller label. This 
correspondence from Morse functions to labelled graphs first appeared on the mathematical 
scene in the work of G. Reeb EH 

A Morse tree has a natural plane structure, i.e. a natural way of linearly ordering the 
branches at every vertex. In particular, every node will have a canonical pair of neighbors 
called successors. Moreover the edges can be canonically colored with three colors called 
North, Northeast and Southeast such every node has a nnique Northern successor, meaning 
that the edge connecting them is a North edge. We called the resulting edge colored plane 
tree a Morse profile because the tubular neighborhoods of two Morse trees in “look the 
same to the naked eye” if they have the same profile. 

In Section 6 we count the number of Morse profiles. We show in Theorem lb.41 that 
their generating function / = f{t) satisfies the cubic equation /(I — tffi = 1 and then we 
determine the Taylor coefficients of this generating function using the Lagrange inversion 
formula. 

In Section 7 we produce a two parameter recurrence for the number of Morse trees (The¬ 
orem EU) which is computationally very effective. In Section 8 we associate to this two 
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parameter family of numbers an exponential type generating function of two independent 
variables and we prove that it satisfies a singular first order quasilinear p.d.e. Fortunately, 
the singularities of this equation can be resolved using a monoidal blowup involving both 
the dependent and the independent variables. The new equation is amenable to the clas¬ 
sical method of characteristics jH] which yields a concrete description of the exponential 
generating function of the number of geometric classes. 

More precisely, if denotes the number of geometric classes of Morse functions with 
2n + 2 critical points and 

_^ ^2n+l 

n>0 ^ ’ 

then in Theorem lO we show that ^ is the compositional inverse of the function 

, r 

, 

\j \— -I- 2st -I- 1 

i.e. ^(0(s)) = s. The proof also reveals a relationship between the Morse function count 
and the classical Weierstrass elliptic functions p. 

The Lagrange inversion formula can in principle be used to compute the Taylor coefficients 
of but the double recurrence in Theorem rm seems easier to implement on a computer. 

The set of topological equivalence classes of Morse functions seems very mysterious at 
this time. It has resisted all our attempts to uncover a computationally friendly structure. 
In section 9 present some partial results. Using some results of Stanley m on counts of 
paths in the Young lattice of partitions we were able to produce a (far from optimal) lower 
bound. There exists at least 1 ■ ■ ■ 3 • • • (2n -|- 1) topological equivalence classes of Morse 
functions on with 2n -|- 2 critical points. We also compute the number of topological 
equivalence classes of Morse functions with at most 10 critical points. 

3. The anatomy of a Morse function on the 2-sphere 

Denote by [n] the set {1,2, ••• ,n}. Suppose / is a excellent Morse function on the 
2-sphere. In this case u(/) is a positive even integer and we set 

Hf) :=^u(/)-l. 

We say that u(/) is the order of /. We denote by Crj C the set of critical points and 
by Aj C K the set of critical values. Suppose u(/) = n, i.e. #Crj = 2n -|- 2. The level 
function associated to / is the bijection 

e : Crf ^ [2n + 2], i{p) ■= #{ q e Crf, f{q)<f{p)}. 

For every i £ [2m + 2] we denote by pi the critical point of / such that i{p) = pi. Note that 
we have a bijection 

i:Af^[2n + 2], Af ^ Cvf [2n + 2], 

We say that a critical value c has level i if £{c) = i. 

For every regular value min/ < c < max/ of / every connected component of the 
sublevel set {/ < c} is a sphere with some disks removed. As we cross a critical value the 
sublevel sets {/ < c} undergoes a change of type Hq, H 2 , or Hf. 

The change Hq corresponds to crossing a local minimum of /. The sublevel set acquires 
a new component diffeomorphic to a closed disk (see Figure |S1 where the new component 
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is depicted in blue). The change H 2 corresponds to crossing a local maximum of / and 
consists of attaching a 2-disk to a boundary component of {/ < c}. 



Figure 3. Attaching a 0-, 2-handle. 

Crossing a saddle point has the effect of attaching a 1-handle. This can be done in only 
two ways: H^, (see Figure 0J. The case corresponds to the gluing of the 1-handle 
to the same boundary component of the sublevel set. The change corresponds to the 
attachment of the 1-handle to different components of the sublevel set. A priori, there is 
another possibility, that of attaching the 1-handle to the same component of {/ < c} but 
in different boundary components. This move can be excluded since it produces a 1-cycle 
which cannot be removed by future handle additions (see Figure El). 



Figure 4. Admissible 1-handle attachments. 
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Figure 5. Inadmissible 1-handle attachment. 


To every regular value min f < c < max / we associate a vector 

b = b{c) = {bo{X'^),h{X^) + l), X<^ = {f<c}. (3.1) 

As we cross a critical value the vector b undergoes one of the changes below 

6^5 + (l,0), 5^5+(0,-l), 




H 


5 + (0,1), b —>5 + (—1,0). 


We denote by m the number of local but not global minima, by M the number of local 
but not global maxima, by n± the number of handles. Note that these integers are 
constrained by the equalities 


The Morse polynomial of / is therefore 


m + n_ + + M 

= 2u 

m — U- — + M 

= 0 

n_ 

= m 

n+ 

= M 

therefore 


(1 + m) + ut + (1 + u 

— m)t^ 


(3.2) 


4. Counting homology equivalenge glasses 

We denote by h{n) the number of ~/i-equivalence classes of Morse functions / : 5*^ —> M 
such that u(/) = n. 

h{z) =Y^h{n)z'^ G Z[[z]]. 

n>0 

We would like to compute h{z). 

Let / G and suppose its critical values are cq < ci < • • • < C2p+i. For i = 0, • • • , 2u 
and 0 < e < min(cj+i — Ci) we set 

di := Ci -h e, 0 < z < 2u. 

The real numbers (di) define a slicing of /. To the function / we now associate using (13.1 II 
the sequence of lattice points 

Pi = b(di) = {bo,i,bi^i + 1) G Z>o, bk^i := dimHk({f < di}, M). 

Set Cl = (1,0), 62 = (0,1). Observe that Pq = P2P = (1,1) and for every i we have 

Pi+i — Pi ^ {±ei, ±62} =: S. 
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We regard the sequence J’(/) = {Po,Pi,--- ,P2 p} as a path in the interior of the first 
quadrant such that every step Pi —> Pj+i has length 1 and is performed in one of the four 
possible lattice directions at Pj, East, North, West, South (see Figureinj. Note that 

/ S' = ‘^g- 


North 





















C 

1 




-► 


East 


Figure 6 . A good path of type consisting of 4 steps East and 4 steps North. 

Define now a lattice path to be a finite sequence of points 7 = {Pq, • • • ,Pl} such that 
Pj+i — Pi € S. The integer L is called the length of the path and it is denoted by L{y). 
The path is called good if all the points Pj are in the interior of the first quadrant. A path 
is called bad if it is not good. For every integer n we denote by T(P, Q; n) the set of paths 
of length n starting at P and ending at Q. Tg and Tfe denote the subsets consisting of good 
and respectively bad paths. We set 

iV(P, Q; n) = #T(P, Q- n), iV,(P, Q; n) = #fP,(P, Q; n), 

W(P,Q;n) = #n(P,Q;n). 

The above discussion shows that 

h{n) = NgiC,C-,2n), C=(l,l). 

The number of good paths is computed in [7j. Consider two points P, Q in the first quadrant. 
Denote by ri, r 2 : the reflections 

ri{ti,t 2 ) = (-DP 2 ), r 2 (ti,X 2 ) = {ti,-t 2 ). 

These reflections generate the Klein group 

K = { l,ri,r 2 ,r = rir 2 } = Z/2 x Z/2. 

Then as in the proof of [71 Thm. 1] we deduce 

W(P, Q- n) + Nb{P, r(Q);n) = Nb{P, ri(Q); n) + W(P, r 2 (Q); n). 

Since r{Q),ri{Q),r 2 {Q) are not in the first quadrant we deduce 

Nb{P,p{Q))=N{P,p{Q)), VpG {r,ri,r2} 

and thus 

Ng{P, Q; n) = N{P, Q; n) - N{P, ri(Q); n) - N(P, r 2 (Q); n) + iV(P, r(g); n). (4.1) 
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Observing that 

N{P,Q-n) = N{^,Q-P-n) 

we deduce 

Ng{C,C;n) = N{0,0;n) - N{0,-2e2;n) - N{0,-2ei;n) +iV(0,-2(ei +e 2 );n). (4.2) 

For t = and x = {xi,X 2 ) € we set P = Consider the step polynomial 

S(t) :=''^P = ti + t2 + G t 2 , ^2 ^]]- 

xgs 

For every x G define 

e,-: C[[ti, t 2 ,^ 2 “']] ^ c, ^ Ay^y ^ A,-. 

yeT? 

Equivalently, if we regard tj as a complex parameter, tj = we have the integral 

formula 

1 “ 


e^A) = ^ A{t)t-^ := ^ 


\tl\ = \t2\ = l 


A{t)t-^d9ide2. 


For every x G Z^ we then have 


so that 


Hence 


A^(0,x;n) = Qs-[KtT), 

Ng{C, C-n) = ji S(t)"(l -tj-tl + tl4) = ^ §{tritl - l)itl - 1). 

-dt = 


h{z^) = / ^ 




— Z§i 


1 — Z ( tl + ^ + t 2 + ^2 ^ ) 


dt. 


=-.K(z,t) 


Observe that 


We write 


and we deduce 

K{z,t) = 


U = ti+t2, V = tit2 
v{v^ — V? + 2v + 1) _ . 2 


= {v^ — + 2 x + 1 ) • 


V — zu{v + 1) 1 — zu{l + 1/x) 

= {v‘^ -u^ + 2v + 1)-Y^ z^u^{l + v-^ Y 


n>0 


n>0 fc=0 


^ ^ - 7/2 + 2x + l)u^v 


n„—k 


^ ^ (fc) + 2u^v^-^ + u^v-'^). 

u—r\ 


n>0 k=0 

Fortunately, very few terms in this sum contribute to (lOl) since 

u^v~^ = 0, Vn / 2k. 


(4.3) 
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Moreover 


We deduce 


E 




j> K{z,t) 


y2m^2-k _ J y2m+2^-k ^ ^ + 6 


This shows that for every m > 0 we have 


h{m) = 


2 m \ ( 2 m\ ( 2 m + 2 \ / 2 m 


m + 2 J \ m 


m + 1 J \m + 1 


2 m\ ( 2 m 


m J \m + 1/ \ m J \ m 


2 m\ f 2 m 


Now observe that 


a + c = 


c + d = 


2 m \ /2m 


m + 2 J \ m 
f 2 m\ ( 2 m 


2 m\ ( 2 m 


m J \m + 1 


2m\ /2m + 1 
m J \ m + 2 


\m J \m + 1/ \m J \ m J 
The sum of the right-hand-sides of the above equalities is 


2m\ /2m\ /2m\ /2m -|- 1 

m J\m J \m J \m + 1 


(X “t“ 2c ~\~ d — 


2m\ /2m + ^2m\ /2m + _ /2m\ /2m -|- 2 

m j \ m + 2 j \ m j \ m + 1 j \ m J \ m + 2 


Hence 


h{m) = a + 2c + d — b = 


m J \ m + 2 


2 m'\/ 2 m + 2 '\ / 2 m + 2 '\/' 2m 

m J \ m + 2 J \ m + 1 J \m + 1 


\m J \ m + 2 J \m + \ J \m + 1J 
1 f 2 m + 2 \f 2 m\_^ 

{m + 2){m + l\m + l )\m J 

where Cm = (^) denotes the m-th Catalan number. We have thus proved the following 

result. 


Theorem 4.1. 


h{m) = Cm+iCm, Vm > 0. 


Example 4.2. Let us test the validity of the above formula for small values of m. We have 
Co = 1, Ci = 1, C 2 = 2, C3 = 5, C 4 = 14, etc. 

We deduce 

h{l) = 2 , h{2) = 10 . 

The two good paths of length 2 from C to C are 

{E,W}, {IV, 5}. 

To list the ten good paths of length 4 it suffices to list only the five good paths whose first 
step points East. The other five are obtained via the reffection in the diagonal xi = X 2 
given by 

E i —> N, W <—^ S. 
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The five good path with initial Eastbound step are 

{E,W,E,W}, {E,E,W,W}, {E,N,S,W], {E,W,N,S}, {E,N,W,S}. □ 

5. Invariants of geometric equivalence classes 

In this section we describe the main combinatorial invariants associated to a geometric 
equivalence class of Morse functions. 

Suppose / is a Morse function / : —> R of order n with critical points {pi, ■ • • ,p2n+2} 

labelled such that 

f{Pi) < f{Pj) i < j- 

In the sequel, by a labelling of a graph will understand an injection from the set of vertices 
to the set of real numbers. For every graph T and every labelling (/? of T we define the level 
function associated to ip 

:= #{ It G V(r); (p{u) < ip{v)}. □ 

We will construct inductively a sequence of labelled forests Fi,--- ,F 2 n+i-,F 2 n +2 which 
replicates the changes in the level sets := {/ < c*, -t- e }, e < minfc(cfc+i — c*,). 

The vertices of F^ are colored in black and white. The black vertices of F^ are bijectively 
labelled with the labels {1,2, ■ • • , k}. The white vertices of Fk have degree one and corre¬ 
spond to the boundary components of the level set Lk which we know is a disjoint union of 
spheres with disks removed. 

The construction of this sequence of forests goes as follows. Fi is the up tree U depicted 
in Figure [7| with the black vertex labelled 1. 


I I A Y 

u D H- H+ 




Figure 7. Simulating handle attachments. 

The passage F^ —> iY+i depends on the topological change Tfc+i, i.e. on the nature 

of the handle attachment as we cross the critical value of level k + 1. As explained in Section 
3, there are four types of handle attachments H 2 , Ff^. 

• i7o-change. F^+i = iY U [7. The new black vertex is labelled k + 1. The labels of the 
previous black vertices are not changed. 
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• i^f2-cliange. Glue the white vertex of the down tree D in Figure [T] to a white vertex of 
corresponding to the component of dL^ where this two handle is attached. The new black 
vertex is labelled k + 1. 

• i7]“-change. Glue the two lower white vertices of the tree H~ in Figure [7| to two white 
vertices of situated in different path components. Label the new black vertex k + 1 and 
ignore the two white vertices where the gluing took place. 

• iL^-change. Glue the lower white vertex of H~ to the white vertex of F^ corresponding 
to the component of dL^ where this 1-handle is attached. Label the new black vertex k + 1 
and ignore the white vertex where the gluing took place. 

The last forest, ^2^+2 consists of single a tree. We will refer to it as the labelled tree 
associated to / and we will denote it by Fj. We denote hy if -.Tf [2n + 2] the associated 
labelling. 

We say that a vertex v is higher than a vertex u (or that u is lower than v) if £/(u) > if{u). 
This tree has 2n -|- 2 vertices, and each of them has either one or three neighbors. We will 
refer to the vertices of degree 3 as nodes. Note that each node has a higher neighbor and a 
lower neighbor. 

To reconstruct the functions from the labelled tree represent it as a tree embedded in 
so that the following hold. 

• The z coordinates of the vertices are equal to the labels. 

• The edges are smoothly embedded arcs. 

• The restriction of the height function z to any edge is an increasing function, or equiva¬ 
lently, the arcs representing the edges have no horizontal tangents. 

We identify with the boundary of a thin tubular neighborhood of the tree in 
(imagine the edges made of thicker and thicker spaghetti). The Morse function is then the 
restriction of the height function z to this surface. 

Definition 5.1. (a) A normalized Morse tree order n is a pair (F,£), where F is a tree with 
vertex set V(r) of cardinality 2n -|- 2 and i : V(r) —>■ [2n + 2] is a labelling such that the 
following hold. 

(al) Every vertex of F has degree one or three. We will refer to the degree 3 vertices as 
nodes. 

(a2) Every node has at least one lower neighbor and at least one higher neighbor. 

A degree 1 vertex is called a maximum/minimum if it is higher/lower than its unique 
neighbor. We denote by N(r) C V(r) the set nodes, and by M^(r) the set of max¬ 
ima/minima. 

(b) Two normalized Morse trees (Fj,.^*), i = 0,1 are said to be isomorphie, (Fo,^o) — 
(ri,.^i), if there exists a bijection 

^:V(ro)^v(ri) 

such that u,v (z V(ro) are neighbors in Fq iff 4>{u),4>{v) are neighbors in Fi and 

ii{(j){u)) = io{u), Vu G Fq. □ 

Observe that the order of a (normalized) Morse tree is precisely the number of nodes. 
The labelled tree {Tf,if) we have associated to a Morse function / on 3“^ is a normalized 
Morse tree and we see that 


/ ~9 5 


(r/,^/) = (r„4)- 


(5.1) 


MORSE FUNCTIONS OF THE TWO SPHERE 


13 


We denote by £(r) the set of edges. Since F is a tree with 2n + 2 vertices we deduce 

#£(r) = 2n + l. 

Definition 5.2. A Morse tree is a pair (r,(/j), where F is a tree and is a labelling of F 
such that (F,£^) is a normalized Morse tree. □ 

Note that the set of Morse trees is equipped with a natural involution 

(F,(/?) <—^ (F,-¥?). 

We will refer to this as the Poincare duality. We will also denote the Poincare dual of a 
Morse tree F by F. A Morse tree will be called self-dual if it is isomorphic to its Poincare 
dual. 

Remark 5.3. (a) The above labelled tree is sometime known as the Reeb graph associated 
to the Morse function /, 0 n cn Cl]. In CD Reeb associates a graph to an excellent 
Morse function / : M —> M, M smooth compact manifold. More precisely one can define 
an equivalence relation on M by declaring two points xq,xi G M equivalent if and only if 
they he in the same path component of some level set f~^{c). The quotient space equipped 
with the quotient topology is then homeomorphic with a finite one dimensional simplicial 
complex whose vertices correspond to the critical points of /. The vertices are then naturally 
labelled by the value of / at the corresponding critical point. When / is an excellent Morse 
function on S'^ then its Reeb graph coincides with the Morse tree defined above. 

(b) It is very easy to produce weaker invariants of a geometric equivalence class. Given 
u,v G M“(F) there exists a unique path in Fj connecting the vertices u,v. mf{u,v) is the 
level of the highest vertex along this path. mf{u,v) is the Mountain Pass level of u and v. 
In Figure 13 we have depicted two Morse functions a,b which are topologically equivalent. 
They are not geometrically equivalent because 

mail , 2 ) = 3, mfc(l,2) = 4. 

Given a Morse tree (F, ip) we a have function 

V(F), 


: P(F) 


^ipiv) = lowest neighbor of v. 




f 

Figure 8. Two Morse functions which are not geometrically equivalent, yet 
topologically equivalent. 

In Figure IS] we have depicted two topologically equivalent Morse functions which are not 
geometrically equivalent because 5/(1) = 3 / 4 = 5g(l). Alternatively, the critical points 
{1,3} of / lie in the same path component of {/ < 3 + e} but the corresponding points of g 
do not he in the same component of {g < 3+e}. Note however that m/(l, 2) = mg(l, 2) = 4. 

□ 
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6. Planar Morse profiles 

The only graphical representations of Morse trees we have depicted so far in Figure 
El E are very special. They are planar, cartesian, i.e. the graphs are embedded in the 
cartesian plane (with coordinates {x,y)), the labelling is given by the height function 
y and the edges are smooths arcs in the plane such that the tangent lines along them are 
never horizontal^. 

Not every Morse tree admits such a nice representation. For example, the Morse tree 
depicted in Figure El does not seem to admit such a representation. 

To understand the origin of the tree in FigureEllet us observe that if we remove the edges 
[6, 7] and [5, 8] then we obtain a collection of Morse trees with no i7]“-handles. These trees 
are determined by the “gradient” flow. We define the gradient of a Morse tree (F, (p) to be 
the map 

V : V(F) ^ V(r) 

where V(i;) = u if u is a local minimum while if u is a node then V (v) is the lowest neighbor 
of V which lies below v. Thus if v is not a local minimum then the edge [u, V(i;)] is the fastest 
descent edge at v. The gradient flow is the (semi-)flow on V(F) defined by the iterates of 
V. Every orbit ends at a local minimum and the set of orbits ending at the same local 
minimum determines a subtree like in Figure El 

However, every Morse tree admits a canonical structure of plane tree. 

12 


2 

1 




Figure 9. A non planar Morse function? 


Proposition 6.1. Every normalized Morse tree (Fjt') of order n has a natural strueture of 
planted trivalent plane tree. i.e. a rooted plane tree such that the root has degree 1, and all 
other vertiees have degrees 1 or 3. with 2n + 2 vertiees. We denote it by 

Proof Indeed, the root corresponds to the global minimum. There is only one edge 
starting at the root. Label this edge by 1. Inductively, once we reached a vertex, then the 
unlabelled edges at this vertex can be ordered by the value of the Morse function at their 
other endpoints. □ 


Q 

This means that the restriction of the height function along any arc has no critical point. 
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As explained in m Exer. 6.19.f] the number of such trees with 2n + 2 vertices is the 
Catalan number . For example the planted ternary plane tree associated to 

the Morse tree in Figure 0 is depicted in Figure Uni 



Figure 10. A plane binary tree associated to a Morse function. 

We want to associate a more refined structure to the Morse tree, called the planar Morse 
profile. This is obtained by drawing the associated plane tree in a very special way. 

Every node v of this plane tree has two successors, a left successor, l{v) and a right 
successor, r{v), where the left successor of v is defined to be the successor with highest 
label. We will draw the edge connecting v to l{v) by a vertical arrow pointing North. The 
edge [u,r(u)] will be represented by an arrow pointing Northeast or Southeast. The arrow 
will point Northeast if the level of r{v) if greater that the level of v, and will point Southeast 
if the level of r(v) is smaller than v. The only edge at the root will be vertical, and will 
point upward. After we do this then we remove the labels indicating the levels. 

Formally, we are coloring the edges withy three colors, North, Northeast and Southeast 
so that each node u has exactly one northern successor v, i.e. the edge [u, v] is colored 
North. We introduce a partial order on the set of vertices by declaring u < v the unique 
path from u to u consists of North or Northeast edges. This order has the property that 
every node has at least one greater neighbor and at least one smaller neighbor. 

Example 6.2. The planar profile of the tree depicted in Figure^] is described in Figure 
El For the reader’s convenience we have not removed the labels so the correspondence with 
Figure EH is more visible. 

The Morse trees depicted in FigureEhave different profiles, yet the corresponding Morse 
functions are topologically equivalent. If Figure El we have depicted the planar Morse 
profiles with 2, 4 and 6 vertices. □ 


Denote by /„ the number of planar Morse profiles with 2n + 2 vertices. Figure El shows 
that 

/o = l, /i = 2, /2 = 7. 


n>0 


Set 
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Figure 11. A planar Morse proEle. 



Figure 12. The planar Morse profiles with at most six vertices. 


We would like to determine f{t). To achieve this we need to investigate the anatomy of a 
Morse profile F. 

First we denote by u(r) the number of vertices and we set 

u(r) := — - — u = 2u + 2. 

We denote hy p = pr its root and by = pjt its only neighbor. The right wing of the 
Morse profile is the maximal sequence of vertices {uq, • • • ,Us}, where 

vo = Py, Vi+i = r{vi), i = 0, •••,s-l. 

The maximality assumption implies that the last vertex Vg on the right wing is a degree 1 
vertex of the tree. 

A vertex vi on the right wing is called a turning point if either Vi is the last point on the 
path, or Vo is the first point on the path and uq > or 0 < i < s and Vi-i <Vi> Vi+i. 

Lemma 6.3 (Monotonicity). If {uo,-- - is the right wing of T and Vi-i,Vi,Vi+i are 
three conseeutive vertiees along this path sueh that Vi-i > Vi then Vi > Uj+i. 


Proof Indeed the edge [vi,ViJ^i] cannot point Northeast since in that case Vi would have 
three higher neighbors Vi-i,l{vi),Vi+i. □ 
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The monotonicity lemma implies that there exists exactly one turning point on the right 
wing, so that it must have one of the three shapes depicted in Figure IT^ In other words, 
the right wing can change direction at most once. 



Figure 13. The possible shapes of the right wing. 

Suppose that we now cut a Morse profile at the vertex pp. We obtain a configuration of 
the type depicted in Figure El where {uq, • • • ,Us} denotes the right wing, Fj denotes the 
left branch of F at Vj , and Vi is the turning point of the right wing. 



Figure 14. Deconstructing a Morse profile. 

The Morse profile is completely determined by the following data. 

• The the length s of the right wing. 

• The location i £ {0,1, • • • , s} of the turning point where the right wing turns^ Southeast. 

• The Morse profiles Tq, • • • , Fs-i satisfying 

u(Fo) H-u(Fs_i) = u(r) - 2 u(ro) H-h u(Fs-i) = u(F) - s. 

Observe that since deg(i;s) = 1 the tree F* is empty. 

Since the turning point can be chosen in (s + l)-different ways we deduce 

n 

(s + l)/no 

s=l noH- \-ns-i=n—s 

Equivalently, for every n > 1 we have 

n 

/„t” = J^(s+l)f( 

s=l noH-|-)is-i=n— s 

it never turns Southeast then we declare the turning point to be the last vertex along the right wing, 
that is i = s. 
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If we sum over n we deduce 


f{t) - fo = ^{s + l)itfy f = 

,^'N.n V J / 


S>1 s>0 

It is more convenient to rewrite the last equality as 

tf = _^_ 

Thus, if we set w = w{t) = tf{t) and Z = r(;(l — wY we deduce Z{w{t)) = t so that w is 
the compositional inverse of Z. Using the Lagrange inversion formula m Thm. 5.4.2] we 
deduce 


w = 


^ n! 


n>0 


d 


'n—l 


dw 


,n—1 I 


W 


Z{w) 


Now observe that 


w 


w 


Z(w) (1 — w)^ \Z{w)J (1 — 

The coefficient of in (1 — is 

(2n)(2n + l)---(3n-2) 

(n — I)! 


and we deduce 

jn-l 

dw^-i l^=o(^ “ = (2n)(2n + 1) • • • (3n - 2) 

so that 

Wn (2n)(2n + 1) • • • (3n — 2) l/3n —2\ 1 /3n + l\ 

n! n! n\n — 1 ) ^ n + l\ n ) 

Thus 

/o=i. /i=i(])=2 , a4Q=7. 

as we discovered before. We have thus proved the following result. 

Theorem 6.4. The number of planar Morse profiles with 2n + 2 vertices is 

1 /3n + A 

n + l\ n )' 

Moreover, the generating series f{t) = Ylin>o fnt^ is a solution of the cubic equation 

f{l-tff = l. □ 
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7. A RECURSIVE CONSTRUCTION OF GEOMETRIC EQUIVALENCE CLASSES 

We would like to produce a two-parameter recursion formula for the number of normalized 
Morse trees of order n. One of the parameters will be the number of nodes and the other 
parameter will be the level of the lowest node. 

If (r, if) is a Morse tree and i is the associated level function we set 

A(r) := v € 3^(r)}, //(h) ;= A(r) - 1. 

Observe that A(r) > 1 since 1 is a global minimum of i so that 1 G M“(/). Hence /u(r) > 1. 
If the order of F is n then the cardinality of iN'(r) is n. Moreover 2n + 2 G E~^{ f) and we 
deduce 

A(r) < n + 2 /i(r) < n + 1. 

In other words, A(r) is the lowest level of a node. In particular we deduce that the vertices 
1, 2, • • • , /u(r)} must be local minima. This shows that if (T, (p) is an un-normalized Morse 
tree then ;u(r) = m if an only if the first m vertices of T are local minima of p and the 
while the vertex of level (m -|- 1) is a saddle point. 

Denote by Fn{m) the number of normalized Morse trees (r,£) of order n such that 
/u(r) > m, that is Morse trees such that the lowest m vertices are local minima. Set 

fn{m) := Fn{m) - Fn{m + 1). 

Note that /n(RT') denotes the number of Morse trees such that the lowest m vertices are 
local minima while the (m + l)-th vertex is a node. Observe that T„(l) is precisely the 
number of normalized Morse trees of order n and 

Fn{m) = ^ fn{k), fn{m) = 0, Vm > n + 1. 

k>m 

Since a Morse tree of order n has exactly n + 2 extrema and one of them must be the highest 
vertex we deduce 


Fn{m) =0, Vm > n + 1 fn{n + 1) = F„(n -G 1). 


Theorem 7.1. 


m 


fnim) = ( 2 ]Fn-i{m - 1 ) 


E X—> Ilf' 


m f m — 1\ f 2n — m + 1 


mi+m2=m—1 \ ni+n2=n—1 

In particular, when m = n + 1 we have 


2 \ mi / \2ni — mi + 1 
n + 1 


Fn{n + 1 ) = fn{n + 1 ) = 


Fni (mi + l)Fn 2 (m 2 + 1) • 


Fn-i{n). 


Proof For every subset C C M we denote by 3^c the set of un-normalized Morse trees 
(F, p) such that 

+ (V(F)) =C. 

We say that C is the discriminant set of p. Note that Tc = 0 if #0 ^ 2Z. We set 
Tc(^ > "i) := { (r,+) G Tc; M(F)>m}, F{C,m) := > m), 

Tc(+= "i) := { (r,+) G Tc; /u(F) = m}, /(C,m) := #Tc(m = 
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Observe that F{C,m) = F{C',m) if #(7 = #C". We set 

3'n ■■= ?'[2n+2], 3‘n{m) := ?'[2n+2](/^ = m). 

If r € then the level m + 1 vertex of F is a node and thus it can only be of two 

types: negative type (—) if m + 1 has two lower neighbors, and positive type (+), if m + 1 
has a unique lower neighbor and thus two higher neighbors. Correspondingly, we obtain a 
partition 

S'nim) = U 

We set F^{m) := so that 

Fn{m) = F+(m) + F~ (m). (7.1) 

We discuss two cases. 

C+. Suppose (r,!") G Denote by v± the two higher neighbors of m + 1 and by u 

the lower neighbor. We set 

k± = k = £{u), k~^ > k~ > m + 1 > k. 

Denote by the graph obtained from T by removing the vertices u, m + 1 and the edges 
at these points. Denote by T^rn the component containing v± and by its vertex set. 
Define f-i- by setting 

V(f ±) = V(r±^) U {r±}, 

where r± has only one neighbor in D^, the vertex v±. Now define (see Figure IT^ 
<F±:V(f±)^R, £±{r±) =m + l, vfi±{w) = £{w), Vu; G V(F)^^^). 


Denote by M± the subset of {1, ■ • • ,m} consisting of vertices which belong to the com¬ 
ponent T^rn- Let m± denote the cardinality of M±. We have 

AF+ U M_ = Mk = {1, • • • , m} \ {£{u)} m+ + m_ = m — 1. 

Clearly M± C V±. If we set 

= V± \ M± 

we deduce that {A-,A^} is a partition of Lm+ 2 ,n = + 2, • • • , 2n + 2}. 

(F±, {p ±) is an unormalized Morse tree such that its first m± +1 vertices are local minima. 
The discriminant set of (p± is 

C± = V± U {m -|- 1} = A± U U {m + 1}, 

so that 

f± G Tc±(m >m± + 1). 

Observe that (F,£) is determined by the integer A; G {1, ■ • • ,m}, the partition 

Mfc = M+UM_, ^M± = m±, 


the partition 


Lm+2,n — A— U A^ 


and the choices of Morse trees (F±, ip ±) with discriminant sets 

C± = M± U A± U {m -t- 1}, 
such that /r(F±) > m± + 1. The order of F± is 


n± = +m±- 1). 


(7.2) 


MORSE FUNCTIONS OF THE TWO SPHERE 


21 




Figure 15. Cutting a Morse tree along the lowest -handle. 

If a± := #A± we deduce a± + m± = 1 mod 2 and 

a± = 2n± — m± + 1, n+ + n_ = n — 1. 

We have thus produced a surjection 

m 

n n “^ck ( fj , > m_ + 1) X 3 ^ Qk{g > m+ + 1)—(7-3) 

fe=l (C^,C^) 

where the ordered pairs {C^,C\) satisfy 

Ci = uA±U{m + 1}, M’Iu M’I = {!,■■■ ,m}\ {k} 

A- U A^ = {m + 2, • • • , 2n + 2}, + 1 mod 2. 

We have a fixed-point-free involution on the left-hand-side of o defined by the bijections 
formally induced by the sign changes -|- <—> — 

3 ^(jk {fi > m_ -I- 1) X {g > m+ -|-1 )—(h > + 1) x (h > "i- + !)• 

The fibers of g> are precisely the orbits of this involution. We deduce 

fn { m ) = ^Yl > m_-h 1) X > m+-h 1). 

k (C^,C^) 

The integer 1 < A; < m can be chosen in m different ways, the partition of 

{1, • • • , m} \ k can be chosen in {^~Y) ways, then we split n — 1 = n_ -|- n+ and we can 
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choose the partition {A-,A^} of {m + 2, • • • , 2n + 2} in ways, a± = 2n± — m± + l, 

and finally, the Morse trees r± G 9^c±("i± + 1) can be chosen in F^^{m± + 1) ways. We 
deduce 


F+{m) 


m 

~2 


E E 

mo+mi=m—1 no+ni=n—1 


m — 1 
mo 


2n — m + 1 
2no — mo + 1 


Fno{mo + l)Fn^{mi + l), (7.4) 


where for typesetting reasons we used the notations mg/i = m±, ng/i = n±. 




Figure 16. Cutting a Morse tree along the lowest -handle. 

C_. Suppose (r,£) G ih~(m). In this case the vertex m + 1 has two lower neighbors 
1 < ki < k 2 < m and a higher neighbor K > m + 1. Remove the vertex k 2 and the unique 
edge at k 2 . The resulting graph T^ is connected (see Figure d. We can now produce 
a Morse tree {T\F) of order n — 1 as follows (see Figure ITbI) . As a graph, F' is obtained 
from Tm by removing the vertex m + 1 and the two edges (m + 1, R), (m+, ki) and then 
connecting ki and K by an edge. Thus 

v(rO =v(r™)\{m + i} c v(r) 

and we set 

/ := ■^lv(r')- 

The the discriminant set of (! is [2n + 2] \ {m +1, ^2 } and F) > m — 1. The normalized 
Morse tree (T, £) is uniquely determined by the integer /ci < /c2 G {1, 2, • • • , m} and by the 








MORSE FUNCTIONS OF THE TWO SPHERE 


23 


un-normalized Morse tree (F', /) of order n — 1, with discriminant set [2n + 2] \ {m + \,k 2 } 
and satisfying > m — 1. We deduce 


m 


Fn ("i) = ( 2 ]Fn-i{m-l). 


(7.5) 


From (Tmi . (Tt:!) and (USD we deduce 


m 


fn{m) = ( 2 ]Fn-iim - 1 ) 


7Tl ^ 

+T E E 


m — 2n — m + 1 


mi 


mi+m2=m—1 ni+n2=n—1 

which is the first equality in Theorem 17.1 1 


2ni — mi + 1 


Fniimi + l)Fn2im2 + 1 ) 


□ 


Example 7.2. For convenience define Fn{m) = Fn{l) for m < 0 so that = 0 for 

m < 0. Let n = 0. Then we have 

/o(l) = 1, him) = 0, Vm 7^ 1 ^ Fo(l) = 1, Fo(m) = 0, Vm > 1. 

We depicted in Figure ITTI the only Morse function on with 2 = 2 x 0 + 2 critical points. 


Figure 17. Morse functions with 2 and 4 critical points. 
Let n = 1. Then Theorem 17.II predicts 

/i(l) = l, /i(2) = Fo(l) = l^Fi(l) = 2, Fi(2) = l. 

The two Morse functions are depicted in Figure ITTI 
For n = 2, m = 1 we have 

1 


A(i) = iE 


j=0 


2j + l 


F,(l)Fi_,(l) = 4Fo(l)Fi(l) = 8. 


Thus we have 8 Morse functions with 6-critical points such that the second critical point is 
not a local minimum. They are depicted in Figure CHI where (5):=selfdual. 

For n = 2, m = 2 we have 

hi2) = Fi(l) + E (. + - j)) 

j=0 V / 

= Fi(l) + Qfi(2)Fo(1) + Qfo(1)Fi(2) = 2 + 3 + 3 = 8 . 

These functions are depicted in Figure [HU where (S'):=selfdual. 

Finally, for n = 2 and m = 3 we have /2(3) = (2)^1 (2) = 3. These three Morse functions 
are depicted in Figure COl Hence 

F2(1) = /2(1) + /2(2) + /2(3) = 8 + 8 + 3 = 19, ^ 2 ( 2 ) = 11, ^ 2 ( 3 ) = 3. 






Figure 18. Morse functions with 6 critical points, and the second is a saddle point. 



Figure 19. Morse functions with 6 critical points, and the first two are minima. 


In Figure lTHlllDll^ the small Latin characters correspond to the profiles depicted in Figure 

ini 

We list below the numbers of Morse profiles, homology classes and geometric classes of 
Morse functions with 2n + 2 critical points, n < 9. The number of geometric equivalence 
classes were computed using a simple Maple procedure based on the above recurrence. For 
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Figure 20. Morse functions with 6 critical points, three minima. 


the reader’s convenience we have included this procedure in the Appendix to this paper. 


n 

2 n + 2 

Morse profiles 

Homology classes 

Geometric classes 

0 

2 

1 

1 

1 

1 

4 

2 

2 

2 

2 

6 

7 

10 

19 

3 

8 

30 

70 

428 

4 

10 

143 

588 

17,746 

5 

12 

728 

5,544 

1,178,792 

6 

14 

3876 

56,628 

114,892,114 

7 

16 

21,318 

613,470 

15,465,685,088 

8 

18 

120,175 

6,952,660 

2,750,970,320,776 

9 

20 

690,690 

81,662,152 

625,218,940,868,432 


(7.6) 


□ 


The computations in the above example give an idea of the complexity of the above 
recurrence and suggest that it can be better organized. We do this in the next section. 

8 . Generating fungtions 
Consider the main recurrence formula 

'm 


fn{m) = ( 2 ]Fn-iim - 1) 


+ - 


m 


E E 


m — 1\ / 2n — m + I 

mi J y2ni — mi + 1 


Fniimi + l)Fn^{m2 + 1 ). 


ni+n2=7i—1 mi+m2=m—1 

Let us first make the change in variables, m = m + 1. We deduce 

/n(m + l)= ^ Fn-i{m) 


m 


2 n — m 


mij V 2ni — mi + 1 


Fniimi + l)Fn2im2 + 1 ). 


Ill V—> V—> 

+Y E E 

ni+n2=n —1 mi+m2=m 

Next, we introduce new functions 

g{m, n) = fnirn + 1), G{m, n) = Fn{m + 1), 0 < m < n. 
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We deduce g{m, n) = G{m, n) — G(m + 1, n). For n > m > 1 we have the equality 

+ V 


G{m, n) — G{m + 1, n) = 


+- 


m 


E E 


m 


2 

2n — m 


niij y2ni — mi + 1 


G{m — 1, re — 1) 

G{mi,ni)G{m2,n2). 


ni-\rn2=n—l mi-\-m2=m 

Now we make the change in variables 

(rre, re) = (x, X + y) (x, y) = (rre, re — rre), H{x,y) := G{m,n). (8-1) 

With these notations the number of geometric equivalence classes of Morse functions with 
2re + 2 critical points is 

Fn{l) = G{0,n) = H{0,n). 

Then we have 

G{m + 1,re) = Ff(x + 1,y — 1), G{m — 1,re — 1) = H{x — l,y). 

For k = 1,2 we make the change in variables in the double sum 

{mk,nk) = (xfc, Xfc + yk) {xk.Vk) = {rnk,nk - mk). 

Then X 2 =reT, — xi=x — xi, y 2 = re — 1 — yi=x + y — 1 — yiso that 

(rrei, rei) + {m 2 , n 2 ) = {m,n - 1) (xi, xi + yi) + (x 2 , X 2 + y 2 ) = (x, x + y - 1) 

^ xi + X2 = X, yi + y2 = y - 1 . 

Now observe that in the double sum we need to have 

2rei — mi + 1 < 2re — rre => 2xi + yi + 1 < 2x + y. 

2 re 2 — m 2 + 1 < 2re — rre 2(x — xi) + (y — 1 — yi) + 1 < 2x + y 0 < 2xi + yi. 

These inequalities are satisfied if and only if 

{xi,yi) e Rx^y-i-.= {{u,v) el?] 0<re<x, 0<v<y-l}. 

For a point (xi,yi) € Rx,y-i we denote by (xi,yi) its reflection in the center of Rx,y-i, i.e. 

(xi,yi) + (xi,yi) = (x,y - 1). 

The recurrence can now be rewritten as 

"x + F 


ff(x,y) - R(x + l,y- 1) = 
^ X + 1 


ff(x - l,y) 


E 


(xi,yi)eRx,y-l 

We now introduce the new function 

1 


X + 2y 
Xij yxi + 2yi + 1 


ff(xi,yi)H(xi,yi). (8.2) 


H(x,y) := 


H{x,y), u{x,y)=x, v{x,y) = x + 2y + 1, 


u{x,y)\v{x,y)V 

H{x,y = —1) = 0, 0! := 1. 

Observe that 

v{x + 1, y - 1) = X + 2y, v{x - 1, y) = x + 2y, 
r;(xi,yi) + r;(xi,yi) = (xi + xi) + 2(yi + yi) + 1 = x + 2y 
We consider two cases. 
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A. X > 0. If we divide both sides of (18.21) by x!(x + 2y)\ we deduce that for x > 0 we have 

{x + 2y + l)H{x, y) - (x + l)H{x + 1, y - 1) 

= —- 1)2/) + —^ H(xi,yi)H{xi,yi) 

{xi,yi)&Rx,y-i 

B. X = 0. If we divide both sides of (lOl by (2y)! we obtain 

(2y + l)iJ(0,y) - Hil,y- 1) = ^ X] yi)^(0,2/ - 1 - yi)- 


yi=o 


Observe that if we let y = 0 in A we deduce 


1 


H{x,0) = -H{x - 1,0) 

so that H{x, 0) = 2“^. Consider the formal power series 

^{s,t) := ^ H{x,y)s'^ty. 

x,y >0 

If we multiply both sides of A and B by and sum over x > 0, y > 1 we deduce 

(x + 2y + l)F(x,y)s'^t^“^ - ^ (x + l)F(x + 1, y - 


E ^-^Hix-l,y)s^ty-^+ E ^ 


E ^iy^i^yi)H{xi,yi)js^ty ^ 

,y — l ' 


Make the change in variables y = y + 1. Then 

E (a: + 2y + 3)i?(x,y + - E {x + l)H{x + l,y)s'^ty 

= E ^^H{x-i,y + i)s^ty+ E ^^[Y.H{xi,yi)H{xuyMs-ty. 


Now make the change in variables x = x + 1 in the third sum. 

E (ic + 2y + 3)i?(x,y + l)s^t^ - E {x + l)H{x + l,y)s'^ty 

= E ^^^—^H{x,y + l)s^^^ty+ E ^^^\'^H{xi,yi)H{xi,yiUs^ty. 


ai>0,y>0 

We obtain 


yRx 


jds{s‘K - s'Kt=o) +2dtn - ds-K = j^ds{s^‘K - s^‘Kt=Q) + ]^ds{s‘K^). (8.3) 


From the equality 


?{t=o = ?f(s,0) =E2- 


x>0 


2 - s 


^5,(sllf) + 25i:if - ds'K = l-ds{s:K^) + 

t 2 It t 2 — s 


we obtain 
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Multiplying both sides by t we obtain 

ds{s:K) + 2tdt^-tdsn=^dsisJ{‘^) + ^ds{s‘^:K) + l, 1>{(S,0) = (8.4) 

The above equality is a first order quasilinear p.d.e. However the initial condition is 
characteristic (see [HI H.l]) and thus the above initial value problem cannot be solved using 
the method of characteristics. To remove the singularities of this equation we blow it up 
via a monoidal change of coordinates 

s = uv, t = v‘^, ^{u,v) = v!K{uv,v‘^). 


Note that 

x,y>Q 

We deduce 


yX ^x+2y+l 

x\{x + 2y + 1 )! 


E « 


6 — a — 1 \ 




a>0,6>l 


a! • h\ 


V = u = st “K = V sJi = = uv 

ds = {dsu)du + {dsv)d^ = = v~^du, 

^1/2 nil 

dt = {dtu)du + {dtv)d^ = -- —du + —dy = -TT^du + —dy = ^{-udy + vdy 

2 2 2v^ 2v 2v^ 

The equation (jH31) can now be rewritten as 

v~^du{u^) + {-udu + n5„)(n“^0 “ vdu{v~^^) = ^du{uv~^^^) + ^du{u^vC) + 1. 

2 2v 

Multiplying both sides of the above equality by v we deduce 

du{uO + {-uvdu + v^dy){v~^i) - vdui = ^du{uf) + '^du{v?C) + v, 


or equivalently 
Dividing by v we obtain 


-vdui + vdyC = ^du{u^^) + '^duiu^O + v. 


-dui + dyi = ^du{uf) + ^du{u^O + 1 - 


This is a first order quasilinear equation with canonical form 

— (1 + n^ + —)duC + dy^ = —+ n^ + 1 . 

The characteristic vector field of this equation is (see D §7.E] or 0 ii.i]) 

H = —(l + n^H—+ (1 + n^ + — )(9.e. 

Consider the curve s —> 7 ( 5 ) described by the initial conditions 

dj 


(8.5) 


u = s, V = 0, ^ = C('S!0)=0 


ds 


= dy 


Along this curve we have H(s,0,0) = —(1 + ^)du + dy + which shows that the initial 
curve is non-characteristic. 
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The characteristic curves of (IH31) are the solutions of the system of o.d.e.-s 


= -(l + < + ^) 


du 
dt 

dv _ 1 

dt ~ ^ 


( 8 . 6 ) 


^ W - 1 + "W? + 


dt I s I 2 ^ 

The graph of ^ is filled by the integral curves of (IHini) with initial points on 7 , i.e. 

m( 0) = s, u(0) = 0, C(0) = 0. (8.7) 

Consider the function 

0 := 2 + ^ + ? = (^ + 0 (^ + l)- 

Observe that the plane curve t 1 —> is a solution of the hamiltonian equation 

u = —d^h 
i = duh 

Since the energy is conserved along the trajectories of a hamiltonian system we deduce 
h{u,^) = const, along the trajectories of (| 8 . 6 B . Thus the solutions of the initial value 
problem (IHSI) + dSZZl) satisfy 

+ u^'^) + u + ^ = s, v = t. 

We interpret the first equality as a quadratic equation in ^ 

,2 


and we solve for ^ 




2 «+(i + y )? + “-« = » 


-(1 + !^) + sj(l + - 2u[u - s) 

u 


Above, the choice of plus sign in the quadratic formula is dictated by the fact that the 
Taylor coefficients of ^ are positive. Thus 


l + tt^+-^ = 'l/(l + ^)^ - 2u{u - s) = \l ^ -+ 2su + 1 




so that 


Set 


du 

dt 


— u^ + 2su + 1 . 


Ps{u) := —— + 2su +1, 9 = 6, 


(n) := [ 

Jo 

9s{u) = Cs-t, Cs = 6s{s) = [ 

Jo 


dr 


vw 


dr 




Then 
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Remark 8.1. We can write u explicitly as a function of s and t using the Weierstrass formula 
[la §20.6, Ex.2] 

3_1 p'{z) + spiz) 


U = ^s{Cs -t), ^s = 




where p{z) = ps{z) is the Weierstrass function with parameters 92 = ^, 93 = j ■ Q 


From the equality 0 = 0s('u = 0) = Cs — t we deduce that the curve u = 0 admits the 
parametrization 

t = Cs = 9sis). 

Now observe that 




-(1 + ^rr + (1 + - Mu - s) 


2s 


u 


((1 + ^) + ■^/(i + \y — 2u{u — s)) ((1 + ^) + \ly + \y — ‘^u{u — s) 


Hence for tt = 0 we have 


X] t) = s, t = Bsis). 


y>o 


( 22 /+ 1 )! 


We have thus proved the following result. 

Theorem 8.2. Denote by the number of geometric equivalence classes of Morse functions 
on 5^ with 2n + 2 critical points and set 

^2n+l 


n>0 


f2n + l)l' 


Then is the compositional inverse of the function 


e{s) = [ 
Jo 


dr 


° ^ + 2 st + 1 


i.e. f,iO{s)) = s. 


□ 


Example 8.3. The Taylor coefficients of ^(t) can be in principle computed from the above 
formula via the Lagrange inversion formula although this procedure is not as effective as 
the recurrence in Theorem o However, we want to test the validity of Theorem 18.21 on 
special cases. 

As in m, for every formal power series / in the variable x we denote by [x'^jf the 
coefficient of x^ in the expansion of /. The Lagrange inversion formula [H Thm.5.4.2] 
implies 

= ( 8 . 8 ) 

We write 

^4 

Ps = l + rs{t), rs{t) = 2st-y + 

Then we have a binomial expansion 


1 


Ps{t) 1/2 = 1 _ lrs{t) + 22 . 2! 


1-3 ,..2 1-3-5^ ^^^3 1-3-5-7 


rs{t) - 


23.3! 


■rsity + 


rs{ty H- 


24.4! 
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Integrating this equality with respect to t € [0, s] we deduce 


Then 


6{s) l-q{s)' \9{s) 


= :q{s) 


= l + 5q{s)+5-6^ + --- . 


We deduce that 




35 15 19 


e{s) 


= 5[s^]q{s) + 15[s^]q{sr = -- + - = 


40 9 24 


5![t^]^ = 19. 


This agrees with the value computed in (ESI). 
Similarly, we have 

l‘'i{ = i[»“i 


e(») 


and we deduce 


7 


0(s) 


= 7[s%{s) + 7-8[s^]^ + 7-8-9[s^] 


,Qis) 


107 


3! 180 


7![t^]^ = 428. 


This too agrees with the value found in ESI)- 


□ 


9. On THE TOPOLOGICAL EQUIVALENCE PROBLEM 

We were not able to find a computationally satisfactory recurrence for the number of 
topological equivalence classes of Morse functions but we could still describe some interesting 
combinatorial structures on this set. 

A regular sublevel set of a Morse function on a sphere is a disjoint union of 2-spheres 
with small open disks removed. 

The topology of such a disjoint union of holed spheres is encoded by a partition tt, i.e. a 
decreasing function 

TT : Z>o ^ Z>o, vr(i) > 7r(i -t- 1), Vi > 0, 

such that 7r(i) = 0 for all i ;s> 0. The length of the partition is the nonnegative integer. 

i{TT) = max{i; 7r(i) > 0} — 1. 

The weight of the partition is the integer 

i>0 

If n = |7r| we say that vr is a partition on n. We denote by T the set of all partitions 
of nonnegative weight and by the set of partitions of weight n. O denotes the unique 
partition of weight 0. 

To a partition vr of weight w and length i there corresponds a disjoint union of i holed 
sphere, 7r(l) holes on the first sphere, 7r(2) holes on the second sphere etc. Observe that the 
weight of the partition is equal to the number of boundary components of the sublevel set. 

We will describe partitions by finite sequences of positive integers, where two sequences 
are to be considered equivalent if one can be obtained from the other by a permutation. For 
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example {4,4,1} = {4,1,4} denotes a union of two spheres with 4-holes and a sphere with 
1 hole. Following a longstanding tradition we will also use exponential notation to indicate 
partitions. Thus 

= My ,hk} . 

mi mfc 

The possible transitions from one sublevel set to the next correspond to four simple opera¬ 
tions on partitions. Given a partition 

hi > h2 ■ ■ ■ > hg > 0 


these operations are 

{hi,-- - ,he} A [hi,-- - ,hg,l}, 

{hi,-- - ,hi,--- ,hg} 1-^ {hi,-- - ,hi + l,--- , hg), 

{hi, - - - ,hi, - - - , hj, - - - hg{ I > {hi, - - - ,hi—i,hi-\-hj 1, hi-i-i, - - - , hj—i, hj-\-i - - - hg{, 

{hi,-- - ,hi,--- ,hg} {hi,-- - ,hi-l,--- , hg}, {hi > 1). 

We will refer to these transitions as moves. Note that the Hq move increases the length and 
the weight by 1, the H{~ move increases the weight by 1 but preserves the length, the 
move decreases the weight and the length by 1, while the H 2 move decreases the weight by 
one but preserves the length. 

We introduce a simplified set of moves. Consider a partition 

hi ^ h2 ^ - - - ^ hg ^ 0. 

The simplified collection of moves is 

• Type U (up) move 

{hi,-- - ,hi,--- ,hg} ^ {hi,-- - ,hi + l,--- , hg}, 

• Type D (down) move 

{hi,-- - ,hi,--- ,hg} ^ {hi, - - - ,hi-l,--- , hg}, hi > 0. 

• Type X-move^ 

{hi, - - - , hi, - - - , hj, - - - hg} I ^ {hi, - - - , hi hj 1, * * • , hj—i, hj-^i, - - - hg}, hi ^ hj ^ 1. 

We can use these moves to produce a directed graph with vertex set T. The procedure 
is very simple. For every move 

T 9 TT ^ tt' G T 

draw an arrow directed from vr to tt'. We denote by IP this graph. Then the number of 
topological equivalence classes of Morse functions with 2n + 2 critical points is equal to the 
number of directed paths of length 2n from the partition 1 to itself. Equivalently, its is 
the number of directed paths of length 2n -|- 2 from O to itself. We denote this number by 
T2n+2- 

We say that a path in T is constructive if it consists only of the [/-moves. A path is 
called simple if it consists only of the U, D moves. Given vr, tt' G IP we write vr ^ tt' if there 
exists a constructive path from tt to tt'. We have the following elementary fact whose proof 
is left to the reader. 

call this the X move since at this moment I don’t know what to make of it. 



MORSE FUNCTIONS OF THE TWO SPHERE 


33 


Lemma 9.1. The following statements are equivalent, 
(a) TT -< Tx'. 

(h) vr C tt' i.e. |7r| < \tx'\ and Tr{i) < 7x'{i), Vi > 0. 


The order relation tt C vr' on T is known as the Young ordering (see [TH Chap.7]). For 
every /x C A we denote by Z„(A,/x) the number of r-chains from ^ to A, i.e. the number 
sequences (/ii, • • • , Hr) satisfying 


5 ■ ■ ■ S C A. 


we have the formula of Kreweras [HI (see also jH §2.4] or [1.11 Exer. 3.63])) 


Zr{n, A) = det 



- hj + 

- j + r ) 




£ = e(X). 


Then if r = jA] — j/rj — 1 then Zr{fJ., A) is precisely the number of constructive paths from fj, 
to A. We denote this number by A). When /r = O this formula simplifies considerably. 
More precisely (7(0, A) is the number of standard Young tableaux (SYT) of shape A (see 
HH Prop. 7.10.3]). More precisely we have the hook-length formula 


(7(A) = (7(0, A) 


lA]! 

n«eA^(«)’ 


where the above product is taken over all the cells u of the Young diagram of A and h{u) is 
the hook length of u, 


h{u) = 1 + if {cells to the right of n} +ff {cells below u}. 



3 

2 

2 

1 


4 

1 


4 

2 


3 

1 


2 



J 



Figure 21. Young diagrams of weight 3 and 4 and 5. 

Example 9.2. In Figure we have depicted all the Young diagrams of weight 3 and 4 
and two Young diagrams of weight 5. The numbers inside the cells the hook lengths. For a 
Young diagram A we set 

h{X) = h{u) 

u£\ 
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and we will refer to it as the hook weight. In the table below we have recorded the hook 
weight of the diagrams depicted in Figure 1211 


A 

h{X) 

3 

6 

2,1 

3 

l3 

6 

4 

24 

3,1 

8 

2 ^ 

12 

2,12 

8 

l4 

24 

2^1 

24 

3,2 

30 


Let us point out that the number of simple paths of length 2n + 2 from O to itself has 
been computed in m Eq. (39)] and it is 

In particular we deduce 

T 2 n +2 > 1 ■ 3 • • • (2n + 1). 

Example 9.3. (a) Since the X-moves do not affect partitions which have only one part 
> 1 we deduce that for n = 2 we have 

Tg = Sg = 15 

so that there are exactly 15 topological equivalence classes of Morse functions on 5^ with 6 
critical points. 

(b) Note that many of the edges of the graph T are edges of the Hasse diagram of the Young 
lattice {T, C). Remove these edges to obtain a new graph Tx- There is an arrow /x —> A 
in this graph if and only A can be reached from /r by one X move, but cannot be reached 
from /i by a ZD-move. For example, (21) can be reached from (22) by a iA-move but cannot 
be reached by an Y-move so they are not connected in Tx- 


•3,2 


2,2 u 3,1 


2 , 2,1 


Figure 22. The graph Tx- 

Let us compute Tg. A path in T of length 8 from O back to O goes only through partitions 
of weight <4. If 7 is such a path which contains an Y-edge then that edge must be the 
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edge (2,2) ^ 3 depicted in Figure HU We deduce that 

Tg = S's + C(2, 2) X (7(3) = 15x7 + 2x1 = 107. 

(c) Denote by Sk{fJ.\X) the number of simple paths in + of length k from fi to X and by 
Tfc(/i|A) the number of all path in IP from /i to A. Upon inspecting Figure!^ we deduce 

Tio = 5io + C{0% 2,1) X (U(0|3,1) + C(0|3)) 

+C(0|3,2) X (C(0|4)+C(0|3)) 

+C(0|2,2) X ( C{0\3, 1) + (7(0|4)) 

= 9 X 105 + 5 X (3 + 1) + 4 X (1 + 1) + 2 X (3 + 1) = 945 + 20 + 8 + 8 = 981. □ 
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Appendix A. Maple implementation 

We include below a simple but far from optimal MAPLE procedure for computing 
the number of Morse functions based on the recurrence satisfied by the numbers H{x,y) 
described in Section 8. 

>Morse:=proc(a::nonnegint, b:: nonnegint) 

>local m, x, y, A;A[0,0]:=1; m:=a+b; 

>for k from 1 to m do 
>for y from 0 to b do 
>x:= k-y; 

>if y=0 then A[x,y] : =1/(2''x) elif x>0 then 

>A[x,y] : = (l/(x+2=i'y+l))=i'( (x+l)*A[x+l,y-l]+ (l/2)=i'(x+l)=i'A[x-l,y] + (l/2)*(x+l) 
^addC add(A[i,j]*A[x-i,y-l-j],j=0..y-1),i=0..x)) else 

> A[x,y]: = (l/(2*y+l))*( (x+1)*A[x+1,y-l]+ (l/2)*(x+l) 

*add( add(A[i,j]*A[x-i,y-l-j],j=0..y-1),i=0..x)) 

> end if; 

> end do; 

> end do; 

> (a!)*(2*b+l)!*A[a,b]; 

> end proc: 

To compute the number H{x,y) in (18.111 use the command 
>Morse(x,y); 

The number of geometric equivalence classes of Morse functions with 2n + 2 vertices is 
obtained using the command 
>Morse(0,n); 

To deal with large n {n > 25) modify the command line 

> (a!)*(2*b+l)!*A[a,b]; 
to 

>A[a,b] ; 

The procedure will then generate the number 
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